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Abstract: We suggest that Diffuse Correlation Spectroscopy (DCS) 
measurements of tissue blood flow primarily probe relative red blood cell 
(RBC) motion, due to the occurrence of multiple sequential scattering 
events within blood vessels. The magnitude of RBC shear-induced diffusion 
is known to correlate with flow velocity, explaining previous reports of lin- 
ear scaling of the DCS "blood flow index" with tissue perfusion despite the 
observed diffusion-like auto-correlation decay. Further, by modeling RBC 
mean square displacement using a formulation that captures the transition 
from ballistic to diffusive motion, we improve the fit to experimental data 
and recover effective diffusion coefficients and velocity de-correlation time 
scales in the range expected from previous blood rheology studies. 
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1. Introduction 

Light scattering methods have been used to probe the motion of suspended particles for the 
last several decades, in either single [?] or multiple scattering regimes [?]. The latter technique, 
known as diffusing wave spectroscopy (DWS) has been extended to heterogeneous multiple- 
scattering media by Boas et al. [?, ?] and has gained acceptance as a method to measure 
perfusion in bulk tissue under the name of Diffuse Correlation Spectroscopy (DCS) [?]. By 
measuring the intensity fluctuations of light diffusely reflected from tissue, DCS can offer a 
measure of microvascular blood flow and has been successfully validated against other blood 
flow measurement techniques, such as arterial spin labeling (ASL) magnetic resonance imag- 
ing (MRI) [?, ?, ?], Doppler Ultrasound [?, ?], Xenon-CT [?] and fluorescent microspheres [?]. 
Given the three dimensional micro-topography of vasculature, red blood cell motion has been 
expected to have the characteristics of ballistic random flow with a uniform spatial velocity 
distribution [?]. Surprisingly though, the good agreement seen in the validation studies cited 
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above requires modeling particle movement as a diffusion random-walk. Two recent studies 
have probed the characteristics of DCS signals more deeply, looking at whether the pulsatile 
nature of blood flow and/or the presence of extravascular tissue shearing contribute to the ob- 
served diffusive scatterer dynamics. Using parallel detection of independent speckles Dietsche 
et al. [?] have measured the auto-correlation of light intensity fluctuations with 26 ms tempo- 
ral resolution, and noted up to 240% variation in the auto-correlation decay time during one 
heart pulsation, as well as a somewhat stronger curvature of the decay curve at the systolic 
maximum flow compared to the diastolic minimum flow. While, this "super-diffusive" decay 
appears to have a slight ballistic flow quality, the data presented by the authors in Fig. 7 of 
Ref. [?] indicates the scatterer motion remains predominantly diffusive throughout the pulsa- 
tion cycle. Further, Ninck et al. [?], using an ex-vivo artificially perfused porcine kidney model 
has shown that, in the absence of blood, the DCS signal carries the signature of extravascular 
tissue shearing. Nevertheless, during pulsatile blood perfusion, the correlation decay curves are 
approximately described by diffusion even as they vary during the pulsation cycle because the 
contribution of extravascular tissue shearing is small. Taking into account the results of these 
studies, Ninck et al. [?] conclude in their discussion that the discrepancy between expected red 
blood cell ballistic flow and diffusion dynamics measured by DCS remains unexplained, but 
perhaps DCS signals might reflect erythrocyte diffusion in the direction perpendicular to flow, 
as has been observed through particle tracking experiments [?], with a magnitude proportional 
to the flow shear rate. 

In this paper we revisit the assumptions made in obtaining blood flow estimates from DCS 
data. In particular we argue that the occurrence of multiple sequential scattering within blood 
vessels would render the DCS measurements sensitive to relative red blood cell motions. As 
noted above, these motions are diffusive in nature, and their magnitude scales nearly linearly 
with blood flow velocity [?], in good agreement with published DCS studies. We also show that 
an effective hydrodynamic diffusion model capturing the transition between early ballistic and 
subsequent diffusive motion results in a modest, but significant improvement in the fit to experi- 
mental data. In addition to an effective diffusion coefficient proportional to blood flow velocity, 
this model also provides a measure of the particle velocity randomization time scale, a poten- 
tially useful tool for blood rheology studies addressing dynamics faster than the millisecond 
range currently accessible using video microscopy methods. 



2. Methods 

2.1. Dynamic light scattering in tissue 

In the context of DWS/DCS, the normalized temporal electric field E(t) auto-correlation func- 
ticmgi(T) = <E(0)E*(t))/(|E(0)| 2 ) is given by [?]: 



gl (T) = £ P(s)exp[- l / 3 kl(Ar 2 (T)) 



ds (1) 



where T is the correlation delay time, P(s) is the normalized probability distribution of detected 
photon path lengths s, ko = 2n/X, with X the wavelength of the illumination laser, (Ar 2 (r)) 
is the mean-squared displacement of scattering particles in time T, and / = 1 /ji s is the photon 
random-walk step length which is equal to the inverse of the scattering coefficient jU s (thus s/l 
gives the average number of scattering events for a pathlength s). The auto-correlation decay in 
biological tissue has been shown to be nearly completely determined by blood flow -related mo- 
tion [?,?]. Thus, given that red blood cells (RBCs) make up over 90% of the corpuscular content 
of the blood, it is reasonable to assume they represent the dominant source of dynamic scat- 
tering events in tissue. In previous DCS investigations, RBC displacement has been modeled 
as either random flow with ^Ar 2 (r)) = V 2 T 2 , where V 2 is the second moment of the velocity 
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distribution, or as Brownian diffusion with (Ar 2 (r)) = 6D/,T, where is the effective Brown- 
ian diffusion coefficient. Unexpectedly, all studies found that the Brownian motion formulation 
leads to a better fit to experimental data than the random flow formulation. Such an observation 
may be explained by noting that the random ballistic flow model assumes successive scattering 
events occur on scattering centers with uncorrected velocity vectors. This is valid if scattering 
events from RBCs are always separated by scattering from extravascular, comparatively static 
tissue, i.e. if the probability of having two or more scattering events in a single blood vessel is 
low. However, this assumption breaks down if the scattering mean free path / is smaller than 
the blood vessel dimension. In fact, studies of blood optical properties indicate the scattering 
length at the common 785-805 nm wavelengths used in DCS instruments is on the order of 12 
jj,m for typical 40% hematocrit [?], while the mean absorption length (l/jU a ) exceeds 3.5 mm. 
Thus photons entering any blood vessel larger than a capillary most likely undergo multiple 
scattering before exiting. Note that 70% of the total blood volume (and an even higher percent- 
age of RBCs because of low capillary hematocrit [?]) is contained in such larger vessels [?]. 
Therefore the majority of intra-vessel scattering events are likely to be part of a sequential scat- 
tering chain. In conjunction with the higher flow velocities in these vessels and the additive 
nature of photon phase change accumulation, we expect sequential within-vessel scattering to 
dictate the photon decorrelation process. Consequently, the random ballistic flow assumption is 
invalidated and we must instead consider the relative motions of RBCs within a vessel. 

Red blood cell dynamics have been the subject of numerous blood rheology studies. As 
also noted by Ninck et al. [?], video microscopy has been used ex vivo to track the motion of 
hemoglobin-depleted ghost RBCs [?] and whole blood [?] in microchannels, as well as in vivo 
in rat venules [?]. It was found that RBCs undergo shear-induced displacements in the bulk 
flow frame of reference that can be characterized by an effective diffusion coefficient D e g on 
the order of 10~ 5 mm 2 /s, much higher than the Brownian diffusion coefficient expected for the 
RBCs in plasma ~5x 10~ 8 mm 2 /s [?]. Most importantly, D e ff appears to scale linearly with the 
shear rate. We postulate that this mechanism underlies the measurement of tissue blood flow 
using Diffuse Correlation Spectroscopy. 

While the simplified Brownian displacement formulation used in DCS literature 
((Ar 2 (r)) =6D£,t) appears to work well, it assumes that the ballistic to random-walk hy- 
drodynamic transition in the RBC diffusion occurs at time scales shorter than those probed by 
DCS measurements. Since there is no data to support this assumption, we remove it by using 
the Langevin formulation for RBC mean squared displacement [?]: 

(Ar 2 (r)) =6D eff (T-T,(l -exp(-T/T c ))) (2) 

where D e ff is the effective diffusion coefficient and % c is the time scale for the randomization of 
velocity vectors associated with RBC scattering events. By Taylor expanding Eq. (??), it can be 
shown that this formulation of the displacement term describes ballistic motion at short delay 
times, and diffusive motion at long delay times. Note that we are referring here to the short time 
scale ballistic motion contained within any diffusive process (including that of erythrocytes 
in the bulk flow frame of reference), and not to the bulk ballistic motion of erythrocytes in 
vasculature as seen in the laboratory frame of reference. 

2.2. Experimental approach and data analysis 

To evaluate the performance of the proposed model, we have applied it to the infant data pre- 
viously reported by our group in Ref. [?]. Briefly, 11 premature infants, between 28 and 34.5 
weeks gestational age were measured several times at weekly intervals, for a total of 66 study 
visits. At each visit, measurements were obtained in up to seven areas of the head using a 
handheld multi-distance probe with source detector distances of 1-2.5 cm. Both frequency do- 
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Fig. 1 . Comparison of data fit errors using the Brownian diffusion (dotted, FVU=0.48%), 
hydrodynamic diffusion (solid, FVU=0.04%) and random flow (dashed, FVU=0.76%) 
mean square displacement models. 



main diffuse reflectance as well as diffuse correlation spectroscopy data was acquired. A fre- 
quency domain instrument emitting 110 MHz modulated light at 8 wavelengths between 659 
and 825 nm was used to estimate absolute optical properties, while the DCS system employed 
a long-coherence length laser operating at 785 nm (CrystaLaser, Reno, NV, USA) and four 
photon-counting avalanche photodiodes (Perkin-Elmer, Quebec, Canada) connected to a dig- 
ital auto-correlator operating at delay times from 200 ns to 1 s (Correlator.com, Bridgewater, 
NJ, USA). Additional experimental details can be found in Ref. [?]. 

To analyze the DCS data we used the semi-infinite medium correlation diffusion formulation 
[?]: 

r ( n t\ ( ex P(-^i) exp(-jQ- 2 ) 

Gi{P,V = -r— 

An \ r\ r 2 

where K 2 = 3jU a ju£ + CHjU^Aq (Ar 2 (r)), a is the probability of scattering from a moving particle 
(RBC), p is the source-detector separation, T\ = (p 2 -hz^) 1 / 2 , r 2 = (p 2 + (zo + lzb) 2 ) 1 ' 2 , zo = 
(jU a + jUj) an d Zb = 1.76/ n[. (assuming the refractive index of tissue to be 1.35). Note that 
we measure experimentally the normalized temporal intensity auto-correlation function, gi(^), 
related to the normalized field auto-correlation function g\ (t) through the Siegert relation (/3 
is a factor dependent on the coherence characteristics of the light source and collection system) 
[?]: 

g 2 {p 1 T) = l+P{gi{p,x)) 2 ;gdp,T) = 




Gl(P.T) 



Gi(p,0) 



(4) 
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3. Results and discussion 



Figure ?? shows a sample intensity auto-correlation curve g2(?) that exemplifies the different 
fits obtained using displacement formulations corresponding to simplified Brownian diffusion, 
hydrodynamic diffusion (Eq. (??)) and random flow, respectively. For this data, it is clear that 
only hydrodynamic diffusion provides a good match to the shape of the auto-correlation de- 
cay. The decay predicted by the simplified Brownian diffusion appears too "slow", while the 
decay predicted by random flow is too "fast", leading to incorrect estimation of the /3 factor 
as well as of the decay rate. To quantify the quality of the fit, we use the statistical metric 
"fraction of variance unexplained" (FVU), defined as the ratio of the model mean squared error 
to the variance of the experimental data. For the data in Fig. ??, the hydrodynamic diffusion 
model has the lowest residuals (FVU=0.04%), followed by Brownian diffusion (FVU=0.48%), 
and random flow (FVU=0.76%). The average values of FVU over the entire set of experimen- 
tal measurements are 0.36% for hydrodynamic diffusion, 0.46% for Brownian diffusion, and 
2.32% for random flow, respectively. As seen in previous studies, the simplified Brownian dif- 
fusion model is found to have significantly lower fit errors compared to random flow. The same 
direct comparison cannot be made between the simplified Brownian and full hydrodynamic dif- 
fusion models because of their different number of parameters. Instead, we perform a statistical 
F-test to determine if the improvement in the hydrodynamic model exceeds the reduction in 
unexplained variance expected from adding an additional fitting parameter (t c ) : 

p _ (SSEm — SSE Def{ ) DoF De{ f 

SSE DeB {DoF DeS - DoF Db ) 

where SSE is the sum of the squares of the residuals and DoF is the number of degrees of 
freedom (number of correlation time bins (137 in our case) minus the number of model pa- 
rameters (2 for Brownian diffusion (D^ and /3), 3 for hydrodynamic diffusion (£> e ff> T c and /3)). 
For hydrodynamic diffusion to be a better model than Brownian diffusion, the corresponding 
F-number must exceed a critical F-value, which for p < 0.05 is 3.91 (as calculated using the 
f inv function in Matlab (Mathworks, Natick,MA)). This is true for 80.1% of the individual 
measurements from our data set, and the whole set F-number is 33.8, corresponding to a highly 
significant p-value of 4 x 10~ 8 (calculated using f cdf in Matlab). 

When Eq. (??) is substituted into Eqs. (??) and (??), a product forms between the prob- 
ability of scattering from a red blood cell a and the diffusion coefficient. This quantity has 
been used as a "blood flow index" in DCS studies, because the value of a generally can- 
not be estimated independently. Figure ?? shows a scatter plot of the aD eB vs. the aDf, 
values obtained from each of our measurements. We observe an approximate relationship of 
<xD e ff = l.Ql aDy + 2.l x 10~ 7 (mm 2 /s). aD eB has a substantially linear relationship with aD/„ 
with a nearly-zero intercept, indicating both parameters can serve as relative blood flow in- 
dices, but aD e ff is expected to provide a more accurate absolute measure. Encouragingly, by 
assuming a=0.1 [?] for brain tissue, our estimated effective diffusion coefficients fall within 
the same range as those determined from video microscopy blood flow studies (10~ 5 — 10~ 4 
mm 2 /s) [?,?,?]. 

With respect to X c , we observed a range between 0.06 and 7.7 /is, with an average of 1.3 /is 
(however values lower than ~ OAjis are not reliable because of the limited time resolution of 
our correlator). A rough estimation of expected T c values may be obtained from hydrodynamic 
diffusion theory. For a rigid spherical particle the velocity decorrelation characteristic time is on 
the order of T v = pa 2 / rj, where a is the particle size, p is the fluid density and 77 is the fluid vis- 
cosity. Assuming a red blood cell diameter of a = 4 fim, rj = 1 .2 cP [?] and the density of water, 
Tv = 13 /is, within an order of magnitude of our measurements. While the deformable nature 
of RBCs and the complexity of blood flow make this comparison less meaningful, our results 



#145333 - $15.00 USD Received 4 Apr 20 1 1 ; revised 13 Jun2011; accepted 17 Jun2011; published 24 Jun 20 1 1 
(C) 201 1 OSA 1 July 201 1 / Vol. 2, No. 7 / BIOMEDICAL OPTICS EXPRESS 2052 




0 L , , , , , 1 

0 0.2 0.4 0.6 0.8 1 1.2 
aD b (mm 2 /s) x10 " 5 

Fig. 2. Scatter plot of aD e ff vs. aDj, across the entire data set. The dotted line indicates the 
diagonal of the plot (ratio=l). 

do indicate the ballistic motion time scale (T < T c ) of the RBC hydrodynamic diffusion process 
is observable in most DCS measurements. Thus the full Eq. (??) should be used to reduce vari- 
ance in the obtained blood flow velocity estimates and to characterize the diffusive transition 
time scale. To further characterize x c we plot in Fig. ??: a) x c vs. flow velocity, represented by 
aD e ff and b) x c vs. inter-RBC distance (expected to be proportional to the inverse cube root 
of the blood hemoglobin concentration HGB~ 1 / 3 ), for all the measurements where R 2 of the 
hydrodynamic fit was greater than 0.999 (giving us confidence in the estimation of x c ). We ob- 
serve a weak but statistically significant decrease in x c with increased blood flow, as well a weak 
decrease with increased inter-particle distance that does not reach a p<0.01 significance level. 
The inverse proportionality between X c and D e ff (and hence blood flow) could be explained as 
an acceleration of the interaction time scale. It is also expected from the short x Taylor ex- 
pansion of the mean square displacement expression (Eq. (??)): (Ar 2 (r)) T = Tr 6(D e f{/x c )x 2 . 
Assuming the short x ballistic displacement takes the form (Ar 2 (r)) = v 2 T 2 , the early ballis- 
tic velocity v is proportional to y / D e ff/' c c- As a first approximation, one could expect v to be 
fairly constant (i.e. dependent on blood viscosity and temperature, but not on the speed of the 
bulk flow), suggesting an inverse-proportional relationship between x c and D e ff. A similar in- 
tuitive explanation is not apparent for the observed decrease in x c with increased inter-particle 
distance. This trend is likely due to complex blood flow mechanisms. Note though that x c is 
known to be affected by particle deformability in colloids [?], thus it may become a useful 
tool to monitor RBC mechanical properties. Such changes can occur due to physiological and 
pathological mechanism, such as the fetal-to-adult hemoglobin replacement, and disease states 
such as sickle-cell anemia. 
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Fig. 3. Dependence of T c on physiological parameters, (a) dependence on the flow velocity, 
assumed to be proportional to aD e ff ; (b) dependence on collision length scale, assumed to 
be inversely proportional to the cubed root of the hemoglobin concentration. 



4. Conclusion 

In conclusion we have provided a mechanistic explanation for the diffusion-like scatterer dy- 
namics observed in diffuse correlation spectroscopy measurements of biological tissue based 
on multiple-scattering within individual blood vessels. We have also shown that experimental 
data is best fit with a hydrodynamic diffusion formulation that includes both ballistic flow at 
short times and diffusive motion at later times. Both the effective diffusion coefficients and the 
transition time scale are within an order of magnitude of expected values. Finally, our results 
imply that, in addition to tissue perfusion measurements, DCS may be a useful tool for blood 
rheology studies. 
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